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Abstract. There is now unequivocal evidence that the jets in FR I radio 
galaxies are initially relativistic, decelerating flows. On the assumption that they 
are axisymmetric and intrinsically symmetrical (a good approximation close to 
the nucleus), we can make models of their geometry, velocity, emissivity and 
field structure whose parameters can be determined by fitting to deep VLA 
observations. Mass entrainment - either from stellar mass loss within the jet 
volume or via a boundary layer at the jet surface - is the most likely cause for 
deceleration. This idea is quantitatively consistent with the velocity field and 
geometry inferred from kinematic modelling and the external gas density and 
pressure profiles derived from X-ray observations. The jets must initially be very 
light, perhaps with an electron-positron composition. 



1. Introduction 

The morphologic al division of ex tragalactic radio sources into two classes intro- 
duced by Fanaro fF Sz Rilev ( 19741 ) has proved to be remarkably robust. Fanaroff 



& Riley had already noted that the edge-brightened (FRII) sources tend to 
have higher radio luminosities than the edge-darkened (FR I) sources, but there 
is also a dependence on stellar luminosity o f the host galaxy which makes the 
division even cleaner ( Ledlow Sz Owen 19961 ). That said, the FRI class is not 



homogeneous: examples of the range of structures are shown in Fig. [TJ Although 
almost all FRI sources show jets on small scales, some appear to be confined to 
the nuclear regions (e.g. 3C 84; Fig. [1^) while others persist almost to the outer 

edges of the source structure (e.g. 3C 296" Fig. [ JP). 

It has long been surmised (e.g. ISimonlll978l : lFanti et al.lll982l ) that the char- 



acteristic structures of FRI sources result from deceleration and disruption of 
their jets, perhaps triggered by interaction with the surrounding intergalactic 
medium. Turning this perception from a vague idea into a quantitative descrip- 
tion of source dynamics proved to be a slow process, but significant progress has 
recently been made using a combination of deep radio imaging, sophisticated jet 
modelling and an understanding of the external e nvironment s of th e sources from 
X-ray observations. These results bear on five of lBlandfordI (|2007l )'s tasks: map 



jet velocity fields; understand the changing composition; measure jet pressures; 
deduce jet confinement mechanisms and infer jet powers and thrusts. 
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Figure 1. Examples of the radio struc tures of FRI sources from the 3CRR 
catalogue (jLeahv. Bridle k Strori]l200(]( l. 

2. FRI jets as decelerating, relativistic flows 




Figure 2. (a) A plot of jet/counter-jet intensity ratio at the brightening 
point , /j/Jni, against n ormalized core power, Pcn for sources from the B2 sam- 
ple (,Laing et al.lll999f ). (b) Depolarization on the counter-jet (DPcountor-jct) 
and jet (DPjet) sides for sources from the same sample with /j//cj > 4 
(jMorganti et al.lll997f) . 



Deceleration of a symmetric, initially relativistic flow provides a very nat- 
ural explanation for the decrease in jet/counter-jet sidedness ratio /j/Zcj with 
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distance from the nucleus observed in FRI jets (e.g. Laing 19931 ). There is now 
direct evi dence for relativ istic motio ns of jet knots on kpc scales in M87 and 
CenA (iBiretta. Zhou fc Owen 19951 : iHarris et al.ll2007l : iHardcastle et al.ll2003l ) 
and many other results are mostly easily explained if FRI let bases are rela- 
tivistic. For example, ij/Icj is correlated with core fraction ( Laing et al. 19991 . 



Fig. [2^) and the counter-jet lobe depolarizes mor e rapidly with increas ing wave- 
length than that associated with the brighter jet (jMorganti et al.lll997l . Fig. [2)3). 




(b) Total intensity 




1 arcsec 



Figure 3. Images of 3C 296 at a resolution of 0.75 arcsec (jLaing et al.ll200^ . 

(a) degree of polarization; (b) total intensity. 



An example of the characteristically asymmetric appearance of FRI jet 
bases is shown in Fig. ^h). As well as this asymmetry, the main and counter- 
jets also show differences in their transverse brightness profiles, in the sense that 
the main jet is more centrally peaked (e.g. Fig. U]). This is a natural manifes- 
tation of a transverse velocity gradient, with higher flow speeds on-axis. There 
is also a systematic difference in polarization structure between the main and 
counter-jets (e.g. Fig. [3^), in the sense that the brighter jet base has a smaller 
degree of polarization on-axis with more prominent regions of longitudinal ap- 
parent magnetic field. If the jets are intrinsically symmetric and relativistic, 
this effect must be due to aberration, which causes the two jets to be observed 
at different angles to the line of sight in the rest frames of their emitting ma- 
terial and therefore to have different observed polarizations. The combination 
of asymmetries in total intensity and linear polarization provides enough infor- 
mation to break the degeneracy between flow velocity and angle to the line of 
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sight, as well as to constrain the three-dimensional structure of the magnetic 
field, allowing us to make detailed models of the jets. 




Figure 4. Jet sidedness for 3C296, derived by dividin g the total-intensi ty 
image from Fig. [3jb) by a copy of itself rotated by 180° (jLaing et al.ll2006[ ). 



3. Models of FRI jets: geometry, speed and field structure 

We have now modelled the flows i n five FRI radio galaxies: 3C 31, B2 0326+39 , 
B2 1553+24, NGC 3 15 and 3C 296 ("La ing Bridl3l2002al : [Canvin Laindl2004l : 



ICanvin et al.ll2005l : lLaing~a l. 2006). We assume that the jets are precisely 



symmetrical, axisymmetric, relativistic flows and derive their geometries, veloc- 
ity fields, emissivity and magnetic-field component ratios by fitting to deep VLA 
observations. We find that a decelerating jet model fits the observed brightness 
and polarization structures and asymmetries well (e.g. Fig. [5|). 

FRI jets on kpc scales can generally be divided into a flaring region (where 
the expansion rate first increases rapidly and then decreases again) and an outer 
region in which the expansion is uniform. We cannot model the jets within a 
kpc or so of the nucleus, where they tend to be faint and poorly resolved by 
the VLA. Where they first brighten, we derive characteristic on-axis velocities 
/3 = f/c f« 0.8. Rapid deceleration then occurs over distances of 1 - 10 kpc 
(always before the jets recollimate) , after which the velocity either stays constant 
or decreases less rapidly. For all sources except 3C 296, the transverse velocity 
profiles are consistent with edge/on-axis ratios ~0.7 everywhere (although we 
cannot rule out evolution from a top-hat profile at the initial brightening). In 
3C 296, the only source with a bridged twin-jet structure that we have observed 
in sufficient detail (Fig. [If), the velocity falls to a low fractional value ^0.1 at 
the jet edge. The jets in 3C 296 may be embedded within the lobes rather than 
propagating in direct contact with the interstellar medium of the host galaxy, 
as appears to be the case for the other sources. Velocity fields for four of the 
sources are shown in Fig. O 

The data hint at further complexity in the jet velocity fields before the 
initial brightening. Measurements of both sidedness ratio and proper inotion 
on pc scales in NGC 315 indicate an apparent acceleration ( Cotton et al. 19991 . 



Fig- CI) • Our VLA observations also show slightly smaller sidedness ratios in the 
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(a) Data 




p = 1 10 arcsec 



(b) Model 
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Figur e 5. A comparison between data and model for 3C296 (jLaing et alj 
|2006[ ). Vectors with magnitudes proportional to the degree of polarization, p 
and with directions along the apparent magnetic field are superimposed on 
grey-scales of total intensity, (a) Data; (b) model. 



faint, kpc-scale jet bases than at or immediately after the brightening point for 
the four sources where we can resolve the former. If we observe a jet stratified 
in velocity with inadequate transverse resolution, the speed inferred for its ap- 
parently brightest emission will be influenced by how changes in its transverse 
emissivity profile interact with the (orientation-dependent) Doppler favouritism, 
and may give a misleading account of the true velocity variation along the jet. 
Both the pc-scale and faint kpc-scale regions are barely resolved in width, so we 
cannot yet tell whether the bulk flow really accelerates. 

On average, the largest single magnetic-field component is toroidal. The 
longitudinal component is significant close to the nucleus but decreases with 
distance and the radial component is always the smallest of the three (e.g. Fig. [5]). 
We also find that the toroidal component is stronger relative to the longitudinal 
component at the edge of the jet. The evolution of the field component ratios 
along the jets before they recollimate is not consistent with flux freezing in a 
laminar flow, which requires a much more rapid transition from longitudinal to 
transverse field than we infer. After recollimation, flux freezing is consistent 
with our results, however. 

Before the jets decelerate and recollimate, the quasi-one-dimensional adi- 
abatic approximation (together with flux freezing) is grossly inconsistent with 
the emissivity evolution inferred from our models. After deceleration and recol- 
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Figure 6. Grey-scales of je t speed deriv ed from relativist ic models for: (a) 
3C31 (iLai ng fc Bridle 2002a): (b) 3C2 96 (iLaing e t al. 20061: (c) B2 0326+39 
(ICanvin fc Laina.2004 ): (d) NGC315 (|Canvin et al..i2005.) . The range is - 
0.8c. 
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Figm'e 7. A comparison of velocity estimates for NGC315 on pc and kpc 
scales, plotted against distance from the nucleus. Filled squares: veloci- 
ties from prope r motions; open sq uares: velocities from jet/counter-jet ra- 
tios (both from lCotton et al1ll999l ): open triangle: velocity derived from the 
jet /counter- jet ratio at 0.4arcsec resolution close to the nucleus. The full and 
dotted lines show m odel fits for the centre and edge of the jets, respectively 
(ICanvin et al.l[200l . 
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Figure 8. Grey-scales showing the fractional magnetic-components in the 
toroidal, long i tudina l and radial directions inferred for the jets of NGC 315 by 
ICanvin at al.l (|2005D . The dominant components are sketched in panel (d). 



limation, however, this approximation describes the emissivity evolution along 
the jets quite well. Details of the radio spectra and high-energy emission are 
given bv lLainget aD (|2007l 1. 



4. Conservation-law analysis and jet energy fluxes 



It is likely that FRI jets decelerate as a result of entraining matter, either 
by injection of mass lost from stars within the jet volume or by ingestion of 
ambient gas from the surrounding IGM via a mixing layer. The kinematic 
models of Section [33 prescribe the variation of velocity along a jet, together 
with its geometry. Given the external pressure and density profiles derived from 
X-ray observations, we can then apply conservation of ene rgy, momentum and 
mass in the quasi-one-dimensional approximation, following iBicknell ( 19941 ) and 
including the effects of buoyancy. Two key assumptions are needed to close the 
problem: the jet is assumed to be in pressure equilibrium with the surrounding 
medium after recoUimation and <I> = He, as expected for an initially relativistic 
jet ($ is the energy flux with the rest-mass contribution subtracted and 11 is the 
momentum flux). In addition, a relativistic equation of state is assumed and 
anisotropic magnetic stresses are neglected. The results of this calculation are 
estimates of the jet energy flux and the variations of internal pressure, density, 
mass flux and entrainment rate as functions of distance from the nucleus. 
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This technique was first apphed to 3C 31 by iLaing fc Bridld (|2002bl ). In 
Fig. [9] we also show the prehminary profiles of pressure, density, Mach number 
and entrainment rate for two other sources - B2 0326+39 and 3C 296 (Laing et 
al., in preparation). The top panels show the internal, external and synchrotron 
minimum pressures. By hypothesis, the first two are matched in the outer 
regions, but we also note that the pressure gradients are quite similar (for 3C 296 
it is possible that the external pressure should be that of the lobe rather than 
the surrounding IGM). The synchrotron minimum pressure is comparable with 
or slightly less than the internal pressure in the outer parts, as expected if the 
latter is dominated by relativistic particles and fields close to equipartition. In 
all three sources, there is a significant overpressure at the brightening point 
(where the analysis starts) and this drives the initial expansion in the flaring 
region. The pressure becomes less than that of the surroundings at the end of 
the flaring region before the jet recollimates and reaches pressure equilibrium 
(the large under-pressure in B2 0326+39 between 2 and 4 kpc is probably due to 
an inaccurate estimate of the outer velocity, which is poorly constrained). The 
jet densities are very low (typically equivalent to 1 proton m~^), and the initial 
density ratios with the surroundings are ~ 10-5. The jets are transonic, with 
Mach numbers ~1. The entrainment rate for 3C31 continues to rise at large 
distances whereas the proflles for B2 0326+39 and 3C 296 both show peaks in 
the flaring region before dropping to very low values. The reason is that the 
kinematic model for 3C 31 requires continued deceleration after recollimation, 
whereas those for the other two sources indicate an approximately constant 
asymptotic velocity. Differences in external density and/or shielding by lobe 
material may be responsible. The entrainment rate profile for 3 C 31 also shows 
a pred iction for the mass input from stars in the jet volume (jLaing &: Bridle 
2002bl ). This is comparable with the required entrainment rate close to the 



start of the model, but not at large distances, suggesting that mass injection 
and ingestion are both important. 

The derived energy fluxes are $ 1.1 X lO^'^, 7 X 10^^ and 1.6 x lO^^W 
respectively, for 3C31, 3C296 and B2 0326+39. It will be o f interest to com- 



pare these estimates with values derived from cavity dynamics (jMcNamara et al 



|2007| . and references therein). Estimates of the mass density and the number 
density of r adiating electrons can be combined to constrain the composition of 
the jets (see Laing Sz Bridleir2002b for a detailed review of assumptions). At the 



brightening point, limiting cases include: 



1. A power-law Lorentz-factor distribution 71(7)^7 oc ^~'^-'^d'y of relativistic 
electrons, each associated with a (cold) proton. In this case, there must 

be a low -GIlGrgy cut-off Sit ^min 

35 (3C 31 and B2 0326+39) or ^ 300 

(3C 296). 



2. An electron-positron jet with entrained thermal matter providing essen- 
tially all of the mass. 

Unfortunately, it is not yet possible to distinguish between these (and other) 
compositions. 
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Figure 9. Profiles of pressure, density, Macfi number and entrainment rate 
along the jets of 3C 31, B 2 0326+39 and 3C 296 f rom a quasi-one-dimensional 
conservation-law analysis (jLaing fc Bridlel[2002U Laing et al., in preparation). 
In the pressure plots, the full, dashed and dotted lines show the internal, 
external and synchrotron minimum pressures, respectively. The dotted curve 
in the entrainment rate plot for 3 C 31 is an estimate of mass input by injection 
from stars within the jet volume (jLaing fc Bridle, 2002b. ). 



5. Where next? 



The techniques reviewed here have led to a detailed picture of the evolution of 
the FRI jets as they flare and decelerate. They work because we can compare 
both jets in the same source: sizes match the resolution of our most sensitive 
instruments, surface brightnesses are (just) high enough, Faraday rotation is 
easy to correct and jet speeds are modest. The key project when EVLA and 
eMERLIN become operational will be an attempt to apply the same methods to 
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other types of relativistic jet: closer to the launching scale and in FRIT sources 
or microquasars. There are also more immediate questions, including: 

1. What makes FRI jets brighten so abruptly in the radio band at distances 
of a few kpc from the nucleus? 

2. Can we decide what combination of stellar mass- loss (injection) and ex- 
ternal entrainment (ingestion) operates in FR I jets? 

3. Can we (finally) detect internal depolarization and hence estimate densities 
of thermal matter now that we understand foreground Faraday screens 
much better? 

4. Do different methods for estimating FRI jet powers (conservation-law, 
cavity dynamics, bow-shock physics) give the same answer? Can they be 
made consistent with radiation-loss timescales? 

5. Do jets really accelerate on pc scales and/or as they fiare or is this effect 
an artefact of inadequate transverse resolution? 
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